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We calculate the thermodynamic functions of Yang-Mills theory to three-loop order using the 
hard-thermal-loop perturbation theory reorganization of finite temperature quantum field theory. 
We show that at three-loop order hard-thermal-loop perturbation theory is compatible with lattice 
results for the pressure, energy density, and entropy down to temperatures T ~ 2 — 3 T c . 
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The goal of ultrarelativistic heavy-ion collision exper- 
iments is to generate energy densities and temperatures 
high enough to create a plasma of quarks and gluons 
called the quark-gluon plasma. One of the chief theoreti- 
cal questions which has emerged in this area is whether it 
is more appropriate to describe the state of matter gener- 
ated during these collisions using weakly-coupled quan- 
tum field theory or a strong-coupling formalism based 
on AdS/CFT correspondence. Early data from the Rel- 
ativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Labs indicated that the state of matter created 
there behaved more like a fluid than a plasma and that 
this "quark-gluon fluid" was strongly coupled 

In the intervening years, however, work on the pertur- 
bative side has shown that observables like jet quench- 
ing [ij and elliptic flow ^ can also be described using a 
perturbative formalism. Since in phenomenological ap- 
plications predictions are complicated by the modeling 
required to describe, for example, initial state effects, 
the space-time evolution of the plasma, and hadroniza- 
tion of the plasma, there are significant theoretical un- 
certainties remaining. Therefore, one is hard put to con- 
clude whether the plasma is strongly or weakly coupled 
based solely on RHIC data. To have a cleaner testing 
ground one can compare theoretical calculations with re- 
sults from lattice quantum chromodynamics (QCD). 

Looking forward to the upcoming heavy-ion experi- 
ments scheduled to take place at the Large Hadron Col- 
lider (LHC) at the European Laboratory for Particle 
Physics (CERN) it is important to know if, at the higher 
temperatures generated, one expects a strongly- coupled 
(liquid) or weakly-coupled (plasma) description to be 
more appropriate. At RHIC, initial temperatures on the 
order of one to two times the QCD critical temperature, 
T c ~ 190 MeV, were generated. At LHC, initial temper- 
atures on the order of 4 — 5 T c are expected. The key 
question is, will the generated matter behave more like a 
plasma of quasiparticles at these higher temperatures? 

In this Letter we discuss the calculation of thermody- 
namic functions of a gas of gluons at phenomenologically 
relevant temperatures. We present results at leading or- 
der (LO), next-to-leading order (NLO), and next-to-next- 



to-leading order (NNLO) and compare with available lat- 
tice data Q for the thermodynamic functions of SU(3) 
Yang-Mills theory. The calculation is based on a reorga- 
nization of the theory around hard-thermal-loop (HTL) 
quasiparticles. Our results indicate that the lattice data 
at temperatures T > 2 — 3 T c are consistent with the 
quasiparticle picture. This is a non-trivial result since, 
in this temperature regime, the QCD coupling constant 
is neither infinitesimally weak nor infinitely strong with 
g s ~ 2, or equivalently a s — g 2 /(47r) ~ 0.3. Therefore, 
we have a crucial test of the quasiparticle picture in the 
intermediate coupling regime. 

The calculation of thermodynamic functions using 
weakly-coupled quantum field theory has a long his- 
tory [5|. The QCD free energy is known up to order 
<7s log(<7 s ); however, the resulting weak-coupling approxi- 
mations do not converge at phenomenologically relevant 
couplings. For example, simply comparing the magni- 
tude of low-order contributions to the QCD free energy 
with three quark flavors (Nf = 3), one finds that the 
contribution is smaller than the <? 2 contribution only for 
g s < 0.9 (a s < 0.07) which corresponds to a temperature 
of T ~ 10 5 GeV - 5 x 10 5 T C . 

The poor convergence of finite-temperature perturba- 
tive expansions of thermodynamic functions stems from 
the fact that at high temperature the classical solution 
is not described by massless gluonic states. Instead one 
must include plasma effects such as the screening of elec- 
tric fields and Landau damping via a self-consistent hard- 
thermal-loop resummation. There are several ways of 
systematically reorganizing the perturbative expansion 
[6| • Here we will present a new NNLO calculation which 
uses the hard-thermal-loop perturbation theory (HTLpt) 
method 04^] and compare with previously obtained LO 
and NLO results. 

The basic idea of the technique is to add and subtract 
an effective mass term from the bare Lagrangian and to 
associate the added piece with the free part of the La- 
grangian and the subtracted piece with the interactions 
loL However, in gauge theories, one cannot simply 
add and subtract a local mass term since this would vio- 
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late gauge invariance. Instead, one adds and subtracts an 
HTL improvement term which modifies the propagators 
and vertices self-consistently so that the reorganization 
is manifestly gauge invariant 12j. 



Formalism : The Lagrangian density for SU(iV c ) Yang- 
Mills theory in Minkowski space is 
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Here the field strength is F^ v = d^A u - d v + 
ig a [Av, A v ], with A^ an element of the SU(iV c ) gauge 
group. The ghost term £ g h depends on the gauge-fixing 
term £ g f. We use MS dimensional regularization with 
a renormalization scale fi and covariant gauge fixing 
£ g f = — (d^A^) 2 /(2£) where £ is the gauge parameter. 
HTLpt is gauge-fixing independent, therefore, all results 
shown below are independent of the gauge-fixing param- 
eter. The independence of the results on the gauge pa- 
rameter was explicitly demonstrated in general covariant 
and Coulomb gauges in Ref. |]|. 

HTLpt is a reorganization of the perturbation series for 
thermal gauge theories. In the case of Yang-Mills theory, 
the Lagrangian density is written as 
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where A£htl collects counterterms necessary to account 
for additional divergences introduced by £htl- The HTL 
improvement term is 
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where = <9 M + igsA^ is the covariant derivative, 



(l,y) is a light-like four-vector, and 



repre- 



sents an average over the directions of y. The term (j^J) 
has the form of the effective Lagrangian that would be 
induced by a rotationally-invariant ensemble of color- 
charged sources in the Eikonal approximation. The free 
parameter mo can be identified with the Dcbyc screen- 
ing mass, but is not assumed to be txid ~ g s T at leading 
order. HTLpt is defined by treating S as a formal ex- 
pansion parameter and expanding in a power series in 8 
around 8 = 0. This generates loops with fully dressed 
propagators and vertices and also automatically gener- 
ates the counterterms necessary to remove the dressing 
as one proceeds to higher loop orders 

If the expansion in 8 could be calculated to all orders, 
the final result would not depend on mo- However, any 
truncation of the expansion in 8 produces results that 
depend on mo- We will first obtain the thermodynamic 
potential fi(T, ot s , rriD, A 4 , 8 = 1) which is a function of 
the mass parameter tud- A prescription is then required 
to determine mo as a function of T and ct s . The canoni- 
cal way to fix the Debye mass in HTLpt is to require the 




FIG. 1: Comparison of LO, NLO, and NNLO predictions for 
the scaled pressure with SU(3) pure-glue lattice data from 
Boyd et al. [j. Shaded bands show the result of varying the 
renormalization scale /i by a factor of two around fj, = 2ttT. 
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FIG. 2: Comparison of LO, NLO, and NNLO predictions for 
the scaled energy density with SU(3) pure-glue lattice data 
from Boyd et al. [J]. Shaded bands show the result of varying 
the renormalization scale fi by a factor of two around fi = 
2ttT. 



thermodynamic potential to satisfy a variational equa- 
tion; however, at NNLO this results in a complex-valued 
Debye mass which causes the thermodynamic functions 
to also become complex. The same problem occurs in 
QED 9] and scalar theories 10] and in the case of gauge 
theories is most likely due to the expansion we perform of 
the resulting integrals in mr>/T (see next section). Due 
to the complexity of the variational Debye mass, here we 
will use a NLO perturbative mass prescription detailed 
below. This prescription guarantees that the Debye mass 
and hence thermodynamic functions are real- valued at all 
temperatures. 

After having fixed the Debye mass as a function of T 
and a s , the free energy, is obtained by evaluating the 
thermodynamic potential at the appropriate value of the 
Debye mass. The pressure, energy density, and entropy 
are then evaluated using standard thermodynamic rela- 
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FIG. 3: Comparison of LO, NLO, and NNLO predictions for 
the scaled entropy with SU(3) pure-glue lattice data from 
Boyd et al. [J] . Shaded bands show the result of varying the 
renormalization scale /j, by a factor of two around /i = 2nT. 
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Thermodynamic Potentials : In this section we 
present the final renormalized thermodynamic potential 
at orders 5° (LO), S 1 (NLO), and S 2 (NNLO). The LO 
and NLO results were first obtained in |8j and we list 
them here for completeness. The thermodynamic poten- 
tials are computed using a dual expansion in g s and m d 
which assumes that at leading order mrj/T is 0(g s ). We 
then only include terms which contribute naively through 
order g^ This dual truncation is not necessary in 
principle; however, in practice it makes the calculation 
tractable. We note that the expansion in mrj/T does 
not spoil the gauge invariance of our final results since 
the gauge parameter dependence cancels prior to the ex- 
pansion in rri£) /T. 

We do not list the renormalization counterterms nec- 
essary but mention that, as in the case of NNLO HTLpt 
QED Q, only systematic vacuum, mass, and coupling 
constant counterterms are necessary to renormalize the 
thermodynamic potential. The thermodynamic poten- 
tials listed below are gauge invariant. Full details of the 
Yang-Mills NNLO calculation will be presented elsewhere 
[l3j : however, we note that the calculation is similar to 
the one presented in Ref. @. 



Leading order 
potential is Q 



The renormalized LO thermodynamic 



LO 



1 — m 7 



deal 



2 
45 
"T 



30 m 3 n 



, A 7 

log h 7 H 

6 2 2 ' 3 



(4) 



where Jideal = — {N 2 — l)ir 2 T 4 /45 is the free energy of 
an ideal gas of non-interacting gluons, 7 is the Euler- 
Mascheroni constant, and we have introduced the dimen- 
sionless parameters rhrj — mrj / (2irT) and fi = /i/(27rT). 



Next-to-leading order : The renormalized NLO thermo- 
dynamic potential is Q 
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Next-to-next-to-leading order : The renormalized NNLO 
thermodynamic potential is 
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where £ is the Riemann £ function. Note that if the lead- 
ing order Debye mass, m 2 ^ w = N c a s /{i-K), is used for 
the Debye mass in ^ we reproduce the known expansion 
of the Yang-Mills free energy up to order a 
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Moss Prescription : The mass parameter itid in HTLpt 
is, in principle, completely arbitrary. To complete a cal- 
culation, it is necessary to specify mo as function of a* 
and T. Unfortunately, similar to NNLO HTLpt QED 
the variational mass prescription gives a complex Debye 
mass. Here we equate the Debye mass used in HTLpt 
with the hard contribution to the Debye mass obtained 
using dimensional reduction i.e. mn = the giving 
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Results : In Figs. H1I31 we show the N c = 3 pressure, 
entropy, and energy density scaled by their respective 
ideal gas limits as a function of T/T c . The results at 
LO, NLO, and NNLO use Eq. for the Debye mass 
in Eqs. (J4j) , (0, and (|6]), respectively. For the running 
coupling we used the three-loop running (l5j ] and varied 
the renormalization scale by two around /i = 2ttT. 
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For the pressure, energy density, and entropy the con- 
vergence of the successive approximations to the Yang- 
Mills thermodynamic functions is improved over naive 
perturbation theory. For example, using the naive per- 
turbativc approach and comparing the full variation in 
both successive truncations and renormalization scale 
variation, one hnds that at T = 3 T c there is variation 
in the pressure of 0.69 < V/Videai < 1-32 Q, whereas us- 
ing HTLpt there is only a variation of 0.74 < V /Pjdeal < 
0.95. Additionally, at NNLO we see that the ^ = 2ttT 
result for the pressure in Fig. Q] coincides with the lattice 
data down to T ~ 3T C and the energy density and en- 
tropy are compatible with lattice data down to T ~ 2 T c . 

However, for all thermodynamic functions we find that 
the NNLO results represent a significant correction to the 
LO and NLO curves. This is unexpected since the LO 
and NLO bands overlap with one another at all tempera- 
tures shown. In addition, in NNLO HTLpt QED [§] such 
a large correction was not observed. For SU(3) Yang- 
Mills, in order for the LO, NLO, and NNLO bands to 
overlap one must go to temperatures T >7T C . One may 
wonder if there is an error in the NNLO thermodynamic 
potential; however, we are confident in this result because 
firstly, we reproduce the known perturbative expansion 
to order a^ 2 in the weak-coupling limit using Eq. ([6]) 
and secondly, there were highly non-trivial cancellations 
of divergences using only systematically predicted coun- 
terterms during the renormalization procedure. We note 
that in scalar theories the corresponding reorganization 
has been pushed to N 3 LO ll| where it has been shown 
that the N 3 LO result is between the NLO and NNLO 
results, indicating an excellent pattern of convergence. 

Conclusions and Outlook : In this Letter we have 
presented a new result for the NNLO thermodynamic 
functions for SU(-/V C ) Yang- Mills theory using the HTLpt 
reorganization. We compared our predictions with lattice 
data for N c = 3 and found that HTLpt is consistent with 
available lattice data down to approximately T ~ 3 T c 
in the case of the pressure and T ~ 2 T c in the case 
of the energy density and entropy. These results are in 
line with expectations since below T ~ 2 — 3 T c a simple 
"electric" quasiparticle approximation breaks down due 
to magnetic/nonperturbative effects. 

We found that at NNLO the variational solution for the 
Debye mass becomes complex and, as a result, we chose 
instead a NLO perturbative mass prescription. The com- 
plexity of the variational solution may be due to the trun- 
cation in rno/T; however, checking this hypothesis will 
require future work. We also found that there was a large 
correction going from NLO to NNLO indicating that per- 
haps the result is not fully converged. Unfortunately, it 
is impossible to say how much a N 3 LO calculation will 
affect things, so again future work is required. 
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